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ABSTRACT 

In this paper, robust transceiver design based on minimum- 
mean-square-error (MMSE) criterion for dual-hop amplify- 
and-forward MIMO relay systems is investigated. The chan- 
nel estimation errors are modeled as Gaussian random vari- 
ables, and then the effect are incorporated into the robust 
transceiver based on the Bayesian framework. An iterative 
algorithm is proposed to jointly design the precoder at the 
source, the forward matrix at the relay and the equalizer at 
the destination, and the joint design problem can be efficiently 
solved by quadratic matrix programming (QMP). 

1. INTRODUCTION 

Recently, amplify-and-forward (AF) MIMO relay systems 
have gained more and more attention from both academic 
and industrial communities, due to its great potential to im- 
prove the wireless channel reliability [JJ, [2J . For practical 
applications, AF MIMO relay systems are to be adopted in 
future communication protocols, such as Winner Project, 
LTE and IMT- Advanced |2J, to enhance the coverage of base 
stations. 

For transceiver design, joint LMMSE transceiver with 
perfect channel state information (CSI) has been investigated 
in 13 and an iterative algorithm has been proposed. Un- 
fortunately, in practice, CSI is generally obtained through 
estimation and perfect CSI is very difficult to achieve. Robust 
transceiver design, which could mitigate such performance 
degradation by taking the channel estimation errors into ac- 
count, is therefore of great importance and highly desirable 
for practical apphcations. 

In this paper, we consider robust linear transceiver design 
for AF MIMO relay systems under imperfect CSI at both the 
relay and destination. The precoder at the source, the for- 
ward matrix at the relay and the equalizer at the destination 
are jointly designed based on minimum-mean-square-error 
(MMSE) criterion. With the channel estimation errors being 
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Fig. 1. Amplify-and-forward MIMO relay diagram. 

modeled as Gaussian random variables, robustness is incor- 
porated into the optimization objective function by taking 
expectation with respect to the channel estimation errors. 
The joint design problem can efficiently solved by quadratic 
matrix programming (QMP) |7|. Simulation results show 
that the proposed robust algorithm performs better than the 
transceiver design without taking channel estimation errors 
into account. 

The following notations are used throughout this paper. 
Boldface lowercase letters denote vectors, while boldface up- 
percase letters denote matrices. The notations Z , Z^ and Z* 
denote the transpose, Hermitian and conjugate of the matrix 
Z, respectively, and Tr(Z) is the trace of the matrix Z. The 
symbol Im denotes the M x M identity matrix, while 0^/ x Af 
denotes the M x N all zero matrix. The symbol E{.} repre- 
sents the expectation operation. The operation vcc(Z) stacks 
the columns of the matrix Z into a single vector The symbol 
(8) denotes the Kronecker product. 

2. SYSTEM MODEL 

In this paper, a dual-hop amplify-and-forward (AF) coopera- 
tive communication system is considered. In the considered 
system, there is one source with Ns antennas, one relay with 
AIr receive antennas and Nb. transmit antennas, and one des- 
tination with AId antennas, as shown in Fig. [T] At the first 
hop, the source transmits data to the relay. The received sig- 
nal, X, at the relay is 



X — HsrPs + ni 



(1) 
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where s is the A^ x 1 data vector transmitted by the source with 
the covariance matrix R^ — E{ss^} = In, P is the precoder 



matrix with a transmit power constraint, Tr(PP^) < P,, 
with Ps is the maximum transmit power at the source. The 
matrix Hgr is the MIMO channel matrix between the source 
and the relay. Symbol rii is the additive Gaussian noise with 
covariance matrix R„j . 

At the relay, the received signal x is multiplied by a for- 
ward matrix F, under a power constraint Tr(FRxF^) < 
P,., where Rx = E{xx^} and Pr is the maximum transmit 
power. Then the resultant signal is transmitted to the destina- 
tion. The received signal y at the destination can be written 
as 

y = HrdFHsrPs + HrdFm + n2 , (2) 

where Urd is the MIMO channel matrix between the relay 
and the destination, and n2 is the additive Gaussian noise vec- 
tor at the second hop with covariance matrix R„2 . In order to 
guarantee the transmitted data s can be recovered at the desti- 
nation, it is assumed that Ns, Mb,, N^, and Md are greater 
than or equal to N. 

It is assumed that both the relay and destination have the 
estimated channel state information (CSI). When channel es- 
timation errors are considered, we have H^^ — H^r + AH^^ 
and Urd = Srd + AUrd, where H^^ and Tlrd are the 
estimated CSI, while AH^r and AHrd are the correspond- 
ing channel estimation errors whose elements are zero mean 

Gaussian random variables. In general, the Mn x Ns matrix 
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AHsr can be written as AH^^ = SI^Hvi/^sr^lr fSl, where 
the elements of the Mn x Ns matrix Hw.sr are independent 
and identically distributed (i.i.d.) Gaussian random variables 
with zero mean and unit variance. The M^ x M^ matrix S^^ 
and Ns x Ns matrix "^^^ are the row and column covariance 
matrices of AH^^, respectively. The matrix AH^^ is said to 
have a matrix-variate complex Gaussian distribution, which 
can be written as AHsr ^ CMMn,Ns (OmrxtVs, Ssr (8) *J^) 
||6l . Similarly, for the estimation error in the second hop, we 
have AHrd ~ CJ\fMD,NR{OMDxNii,'^rd «) *Jd), where the 
Mjj X Md matrix S^d and Nr x Nr matrix \I/J^ are the 
row and column covariance matrices of AHrd, respectively. 
It is assumed that H^r and Urd are estimated independently, 
so the channel estimation errors, AH^r- and AHrd, are inde- 
pendent. 

3. PROBLEM FORMULATION 



MSE expression (O can be written as 

MSE(G,F,P) 

= EAH,,.AH„{Tr((H,.dFH,rP)(H,dFH,,P)«)} 
+ EAH,,{Tr((GH,.dF)R„,(GHrdF)")} + Tr(GR„,G«) 
+ Tr(Iiv) - Tr(GH,dFH,rP) - Tr((GH,dFH,rP)") (4) 

Because AHs^ and AH^d are independent, the first term of 

MSE is 



EAH,,„AH,,{Tr((GH,.dFH,rP)(GH.dFH,rP)'')} 

= Tr(GEAH,, {H.dFEAH„{H.rPP''H?jF«H?d} G^^). 

(5) 

For the inner expectation, due to the fact that the distribution 
of Atlsr is matrix-variate complex Gaussian with zero mean, 
the following equation holds ||6| 

EAH,,^{H,rPP"Hfr} = Tr(PpH*sr)S,r + Hs.PP^flf, 
= np. (6) 

Applying ^ and the corresponding result for AH^d to (|5]), 
the first term of MSE in (HI becomes 

Tr (gEah,,, {H,dFEAH„{H.rPP"Hf4F'^Hj:'d} G" 
= Tr(G(Tr(FnpF'^*rd)S,d + HrdFRpF^fiJ^JG"). (7) 

With similar calculations applied to the second term of MSE, 
the total MSE in (|4]l can be shown to be 



MSE(G,F,P) =Tr G(H,.dFRxF"Hj:'d + K)G" + Tr(I 



Tr (GH,dFH,,P) 



(8) 



where R.^ 



Ht 



R„, and K = Tr(FRxFH*rrf)Srd 



R„2 . Notice that the matrix Rx is the autocorrelation ma- 
trix of the receive signal x at the relay. Finally, the joint 
transceiver design can be formulated as the following opti- 
mization problem 

min MSE(G,F,P) 



s.t. Tr(PP") < Ps, Tr(FRxF'^) < Pr 



(9) 



At the destination, a linear equalizer G is adopted to detect 
the transmitted data s. The problem is how to design the linear 
precoder matrix P at the source, the linear forward matrix F 
at the relay and the linear equalizer G at the destination to 
minimize the mean square errors (MSE) of the received data 
at the destination: 

MSE(G,F,P) 
= E{||(GH,.dFH,.P - Iiv)s + GH.dFni + Gnall'} (3) 

where the expectation is taken with respect to s, AH^,., 
AHrd, ni and n2. Since s, rii and n2 are independent, the 



4. THE PROPOSED SOLUTION 

In this section, we derive an iterative algorithm to solve for 
P, F and G. In the following, it is shown that given any two 
variables of P, F and G, the remaining one can be efficiently 
solved. Therefore, the proposed algorithm computes P, F 
and G iteratively, starting with initial values. 

Design of G : When the precoder P at the source and 
the forward matrix F at the relay are fixed, the optimiza- 
tion problem (|9]l is an unconstrained convex optimization 



problem for G. The optimal equalizer G must satisfy 
(9MSE(G, F, P)/5G* = which gives 

G = (H,dFH,,P)"(H,dFRxFHH^d + K)-\ (10) 



Design of F : When P and G are fixed, the optimization 
problem (|9]l become s 

minp MSE(G,F,P) s.t. Tr(FRxF") < P^. (11) 

Generally speaking, the optimization problem (fTTT l is a 
quadratic matrix programming (QMP) problem with the vari- 
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able FRx and only one constraint. We can formulate it 
into a semi-definite programming (SDP) problem to solve for 
F. However, because there is only one constraint, in the fol- 
lowing, we introduce another algorithm to compute F based 
on Karush-Kuhn-Tucker (KKT) conditions, and has a much 
lower complexity. The corresponding KKT conditions of ( fTTT l 
are given as follows 



F = (H?dG"GH,d + *,dTr(GS,dG") + XIn^)'^ 

i-H /-iH-nH-fVH 



X H^rfG P H.^Rx . 



A(Tr(FRxF'') - P,.) = 
A > 0, TrfFRxF") < Pr 



(12a) 
(12b) 
(12c) 



Obviously from (I12al i. in order to compute the optimal F, 
the Lagrangian multiplier A should be calculated first. How- 
ever, there is no closed-form solution of A simultaneously sat- 
isfying ( |12bt and ( |12ct . Below we propose a low complexity 
method to solve ( I12bl l and ( I12cb . First, notice that in order 
to have (fT2b] l satisfied, either A = or Tr(FRxF") = P,. 
must hold. If A == also makes jllcj satisfied, A = is a 
solution to (|12bb and (|12ct . Since given G and P, the op- 
timization problem ( fTTT l is a convex quadratic programming 
problem of F, which has only one solution for F, A = is the 
only solution to ( I12bl ) and ( I12cl ) in this case. 

On other hand, if A = does not make ( |12c| ) satisfied, we 
have to solve Tr(FRxF^) = Pr- It can be proved that when 
G and P are fixed, the function /(A) = Tr(FRxFH) is a 
decreasing function of A which satisfies 



< A < \/Tr(HH^GHpHHH R-iH«.PGH,d)/P.. (13) 

Due to space limitation, the proof is not presented here. 
Based on this result, A can be efficiently computed by a 
one-dimension search, such as bisection search or golden 
search. Since Tr(FRxF^) = P^ is a stronger condition than 
Tr(FRxF^) < Pr, ( I12cl ) is satisfied automatically in this 
case. In summary, we take A = 0, if /(O) < Pr, and solve 
/(A) — Pr otherwise. 

Design of P : When F and G are fixed, after a lengthy 
and tedious derivation, it can be shown that the optimization 
problem ^ is equivalent to the following QMP problem |j4l 

mm Tr(P"AoP) + 27^(Tr(BJfP)) + co 

s.t. Tr(P"AiP) + 27^(Tr(Bf P)) + Ci < 

Tr(P"A2P) + 27^(Tr(Bf P)) + C2 < 0, (14) 



where the parameters are defined as follows 

Ao = *,.Tr(FS,,F"M) + H?,F"MFH»„ 

M ^ *,dTr(GE,dG'^ 



H^^G GHrtj, 



Bn 



(GHrdFH 



, Co 



: Tr(G(Ri + R„,)G") + Tr(Iiv), 



Ri = Tr(FRniF 'I'rd)Srd + HrdFR„jF H,,d, 
Ai — Ins, Bi — Ons,n, C\ = — Ps, 
A2 = *«.Tr(FS«.FH) + Hf,F"FH,„ 



B2 



0]vs,iv, C2 = Tr(FR„iF") - P„ 



(15) 



It is known that QMP problems can be transformed into 
semi-definite programming (SDP) problems which can be ef- 
ficiently solved by interior point polynomial algorithms Q . 
Based on the properties of Kronecker product and the follow- 
ing definition 



a = 



I]v®Ai vec(Bi) 
vec^(Bi) d 



i = 0,l,2, 



(16) 



the optimization problem ( fT4b is equivalent to 

min Tr(noX) 

X 

s.t. Tr(niX) < 0, Tr(n2X) < 
X= [vec'^(P) l]^[vec"(P) 1] 



(17) 



If the constraint Rank(X) = 1 is relaxed (it is a well-known 
semi-definite relaxation (SDR) Q, |[8|), we have the follow- 
ing SDP relaxation problem 

min Tr(noZ) 



s.t. Tr(niZ) < 0, Tr(n2Z) < 
[Z]]vjVj,+i,]viVj,+i = 1, Z >; 0, 



(18) 



where Z is a Hermitian matrix. Because the QMP problem 
( fT4] i is a convex quadratic programming problem, the relax- 
ation gap of SDR is zero. In other words, the optimization 
problems ( fTTI i and ( fTSI ) have the same optimal solution H, 

Summary and Convergence Analysis: Initiahze P and 
F which satisfy Tr(PpH) = P^ and Tr(FRxF") = P^. For 
simplicity, we can take P oc I and F oc I. Then the proposed 
iterative algorithm proceeds between ( fTOf ). ( I12ab and ( fTSl l. un- 
til llMSEi - MSEi_i|| < t where MSE^ is the MSB % in 
the zth iteration, and i is a threshold. Since for any two of the 
P, F and G fixed, the optimization problem (|9l) is a convex 
problem for the remaining variable, the proposed algorithm 
is an alternative projection algorithm which is guaranteed to 
converges. 

5. SIMULATION RESULTS AND DISCUSSIONS 

In this section, we will investigate the performance of the pro- 
posed algorithm and for the purpose of comparison, the algo- 
rithm based on the estimated channel only (without taking 



the channel errors into account) O is also simulated. In order 
to solve the SDP problem, the matlab toolbox CVX is used 
IfTOl . In the following, we consider an AF MIMO relay sys- 
tem where the source, relay and destination are equipped with 
samenumberof antennas, i.e., Ns = Mn = Nr = Mo = 4. 
The estimated channels Yisr and H^d are randomly generated 
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Here the channel estimation algorithm in fSl is adopted, the 
correlation matrices of channel estimation errors are in the 
form 13. 

*sr = R.T,sr, Ssr = (Je (Ijv/n + <^e'^R,sr) ) 
*rd = R-T,rd, 5]rd = CTg (Ia/j^ + 0"e R-fl,rd) i 

where o-g denotes the estimation error variance fS]. The 
matrices Rr.sr and R_r, sr are the transmit and receive cor- 
relation matrices in the first hop, respectively, and similar 
definitions apply to Rr.rd and R/j^rd for the second hop. The 
widely used exponential model [5] is chosen for the transmit 
and receive channel correlation matrices, i.e., (R.T,rd)i,j = 

{R.T,sr)i.,j = a)^~-'\ (Rfl,rrf)ij- = (Rfl,sr)i,j = j3\^~^\ 

where a and /3 are the correlation coefficients. 

We define the signal-to-noise ratio for the source-relay 
link (SNRsr) as E^/Ni ^ Ps/Tr(R„^), and is fixed as 
Eg/Ni — 30dB. At the source, four independent data streams 
are transmitted. For each data stream, 10^ independent QPSK 
symbols are transmitted. The SNR for the relay-destination 
fink (SNR^d) is defined as E^/Na = Pr/Tr(R„J. Each 
point in the following figure is an average of 1000 indepen- 
dent realization of estimation errors. 

Fig. |2] shows the bit-error-rate (BER) performance of the 
proposed algorithm and the algorithm based on estimated 
channels only with different a1, when a = 0.5 and /3 = 0.4. 
It can be seen that when the channel estimation errors de- 
creases, the performances of both algorithms improve and 
they coincide at (7^ = 0. Furthermore, the performance of the 
proposed algorithm is always better than that of the algorithm 
based on estimated channels only. 

6. CONCLUSIONS 

In this paper, based on the Bayesian framework, robust linear 
transceiver design for dual-hop AF MIMO relay systems has 
been considered. The precoder matrix at the source, the lin- 
ear forward matrix at the relay and the linear equalizer at the 
destination have been jointly designed based on minimum- 
mean-square-error (MMSE) criterion. An iterative algorithm 
is proposed, and at each step, the design problem can be for- 
mulated as a QMP problem which can be efficiently solved. 
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Fig. 2. The BERs for the proposed iterative algorithm and the 
algorithm based on estimated channels only for different ct^, 
when a = 0.5 and (3 — 0.4. 

Simulation results showed that the performance of the pro- 
posed robust algorithm is always better than that of the algo- 
rithm based on estimated channels only. 
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